Vector current conservation and neutrino emission 
from singlet-paired baryons in neutron stars. 
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Abstract 

Neutrino emission caused by singlet Cooper pairing of baryons in neu- 
tron stars is recalculated by accurately taking into account for conserva- 
tion of the vector weak currents. The neutrino emissivity via the vector 
weak currents is found to be several orders of magnitude smaller than 
that obtained before by different authors. This makes unimportant the 
neutrino radiation from singlet pairing of protons or hyperons. 

One of the mechanisms leading the neutron star cooling, specially for tempera- 
tures near the critical one Tc, consists on the recombination of thermally excited 
baryon BCS pairs into the condensate. This process has received the attention 
of many authors, and is currently thought to be dominant, for some ranges of the 
temperature and/or matter density (modulo the uncertainties arising from the 
incomplete knowledge of the gap value). A better understanding of this process 
is, therefore, of great importance for the secular evolution of such objects. 

Under the description of nuclear matter in the nonrelativistic limit, the sim- 
plest case for baryon pairing corresponds to two particles bounded in the ^5*0 
state. The neutrino emission for recombination into this state was first cal- 
culated by Flowers et al. and reproduced by other authors [2], 0. The 
neutrino energy losses due to pairing of hyperons [5] are also discussed in 
the literature as possible cooling mechanisms for superdense baryonic matter in 
neutron stars. Nowadays, these ideas are widely accepted and used in numerical 
simulations of neutron star evolution [S], [7], |H1- 

In the case of singlet pairing, the averaged weak axial current vanishes, and 
the emissivity is directly mediated by the weak vector current. As it is well 
known, the vector current possesses the property of being conserved by nuclear 
and electroweak interactions. Of course, this fundamental property has to be 



accounted for in any calculation of neutrino emission processes. As we show 
below, however, previous calculations did not pay attention to this particular 
topic. This translates into a dramatic overestimation of the energy production 
from the process under consideration. 

Let us recall shortly the main steps in the above calculations. The low- 
energy Lagrangian of the weak interaction may be described by a point-like 
current-current approach. For interactions mediated by neutral weak currents, 
it can be written as^ ^ 



Here Gf is the Fermi coupling constant, and the neutrino weak current is given 
by = D^^ (1 — 75) v. The vacuum weak current of the baryon is of the form 
Jn — "0 (Cvln ~ CaI^Iz) "0; where ij^ represents the baryon field, and the weak 
vertex includes the vector and axial-vector terms with the corresponding cou- 
pling constants Cy and Ca- 

As mentioned above, in the case of singlet pairing only the vector current con- 
tributes. The nonrelativistic limit for this current is taken as tps^^ipB — ^ ^b^b, 
all others being zero. Here is the second-quantized nonrelativistic spinor 
wave function. The process is kinematically allowed due to the existence of a 
superfluid energy gap A, which admits the transition with time-like momentum 
transfer K = (w, k), with lu = lui + UJ2 and k — ki + k2 being the energy and 
momentum carried out by the freely escaping neutrino pair. We have w > 2A 
and uj > k . 

The relevant input for this calculation is the recombination matrix element 
between the baryon state, which has a pair of quasi-particle excitations of 
momentum-spin labels (p, up; p', down), and the same state but with these ex- 
citations restored to the condensate. To the leading (zero) order in A: <C Pf, 
this matrix element is usually estimated as ^ 

\Mb\' = ^ (2) 

where ep is the quasi-particle energy, as given by Eq. I^. As a result, the 
neutrino energy losses at temperature T < Tc arc found to be: 



where M* is the effective nucleon mass, y ^ /S./T ^ z = \/ + y"^, and Afi, = 3 
is the number of neutrino flavors. 

The naive estimate Q is inconsistent with the hypothesis of conservation of 
the vector current in weak interactions. Indeed, a longitudinal vector current of 
quasi-particles consisting only on a temporal component can not satisfy the con- 
tinuity equation. It is well known, however, that the Bardeen-Cooper-Schrieffer 

^In what follows we use the Standard Model of weak interactions, the system of units h = 
c = 1 and the Boltzmann constant kg = 1. The fine-structure constant is a = e^/47r = 1/137. 
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theory of superconductivity is gauge invariant |H] and that the current conserva- 
tion can be restored if the interaction among quasi-particles is incorporated in 
the couphng vertex to the same degree of approximation as the self-energy effect 
is included in the quasi-particle JOli In the present paper we recalculate 
the neutrino energy losses with allowance for conservation of the weak vector 
current. 

It is convenient to use the Nambu-Gorkov formalism, where the quasi- 
particle fields are represented by two-component objects 

np)-( t^'D- (4) 



Here ipi {p) is the the quasi-particle component of the excitation with momen- 
tum p and spin cr, and ip^ {~p) is the hole component of the same excitation, 
which can be interpreted as the absence of a particle with momentum — p and 
spin —a. The two-component fields obey the standard fermion commutation 
relations 

With the aid of the 2x2 Pauli matrices Ti — (fi, f2, fa) operating in the particle- 
hole space, the Hamiltonian of the system of quasi-particles can be recast as [TT] 

H = Ho + Hi, 

where 

P 

is the BCS reduced Hamiltonian, the nonrelativistic energy is measured rela- 
tively to the Fermi level 

and fi is the Fermi energy. The BCS reduced Hamiltonian bears a resem- 
blance to the one describing the Dirac equation. It has eigenvalues po = ±ep 
with 



a + (6) 

which correspond to excited states in the particle-hole picture, while the ground 
state (vacuum) is the state where all negative energy "quasi-particles" (e < 0) 
are occupied and no positive energy particles exist. The positive and negative 
states are separated by an energy gap 2A. 

The Hamiltonian of residual interaction among quasi-particles has the fol- 
lowing form 
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As follows from the Hamiltonian the inverse of the quasi-particle propagator 
can be written as |11|: 

= Po-epT3-Afi, (7) 

The self-energies are, in general, complex numbers due to the instability of 
single particles. However, to the extent that the single-particle picture makes 
some physical sense, we will ignore the small imaginary part of the self-energies, 
and describe the quasi-particles with the aid of wave- functions. The states of 
quasi-particles obey the equation 

g;-i*p = 0. (8) 

The solution to this equation corresponding to the energy pq — Sp and spin 
state Xa has the following form 





with 



There is also a solution of negative frequency po — —^p 

V "px^ / 

which corresponds to the same excitation energy. This solution is connected to 
the hole state by the particle-antiparticle conjugation 

C : ^ C*^ = f2*^ 

which changes quasi-particles of energy- momentum {pq, p) into holes of energy- 
momentum (— po,— p)- 

The components of the bare vertex 

^ \ if ^ = z = 1,2,3 ■ 

are 2x2 matrices in the Nambu-Gorkov space. As already mentioned, the longi- 
tudinal current corresponding to the bare vertex does not satisfy the continuity 
equation. To restore the current conservation, one must consider the modifica- 
tion of the vertex 7^ to the same order as the modification of the propagator 
is done. The relation between the modified vertex F'' and the quasi-particle 
propagator (O is given by the Ward identity 

K^r^ {p',p) = faQ-i {p) - {p') f3, (12) 
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where K — (i^,k) is the transferred momentum. The plane wave solutions 



Mp,a exp (zpr — lEpt) , Up'^Q' exp {—ipr + iept) 

for 5* and obey the equations Q^^ (p) '«p,a — 0, and Up/^a'Q'^^ ip') — ^■ 
Therefore, the Ward identity implies conservation of the vector current on the 
energy shell of the quasi-particles. Following the prescriptions of quantum elec- 
trodynamics, an approximation which satisfies the Ward identity (and hence 
the continuity equation) is the sum of ladder diagrams. 

Consider first the case of electrically neutral baryons. Then the corrected 
vertex can be found from the following Dyson equation 



r^ip-K,p) = r{p-K,p) (13) 

+' I 7^ iv' - K) [p' - K,p') Q ip') f,V„, , 

J (ZTT) 



(27r 

where the "dressed" particles interact with the same primary interaction Vpp' 
which produces the self-energy of the quasi-particle. 
In the limit K ~ {uj, 0), the Ward identity gives^ 

ro (p-X,p)=f3--if2A 

The poles of the vertex function correspond to collective eigen-modes of the 
system. Therefore, the pole which appears at uj ^ 0, k — implies the existence 
of a collective mode, which plays an important role in the conservation of the 
vector current. The corresponding nonperturbative solution to Eq. H13|l has 
been found by Nambu Jl] (see also ^3): 

To{p-K,p)^f^~2^T2A-^^-^ (14) 

^=^-2^f,A /\ (15) 

The poles in this vertex correspond to the collective motion of the condensate, 
with the dispersion relation lo = ak, where — Vp/3. 

The effective vertex satisfies the Ward identity H12|l , and thus the continuity 
equation on the energy shell 

ujTo - kr ~ (16) 

We are now in a position to evaluate the matrix element of the vector weak 
current. In the particle-hole picture, the creation and recombination of two 
quasi-particles is described by the off-diagonal matrix elements of the Hamilto- 
nian, which corresponds to quasi-particle transitions into a hole (and a corre- 
lated pair). Thus, we calculate the matrix element of the current between the 

■^To obtain the weak vector current this vertex should be multipUed by the weak coupUng 
constant Cy. 



5 



initial (positive-frequency) state of a quasi-particle with momentum p and the 
final (negative-frequency) state with the same momentum p. 



Let us consider separately the contributions from the bare vertex, given by 
the first term in Eq. H14|l . and the collective part, given by the second term, so 
that Mu= + 

Making use of the wave functions described by Eqs. (|9J), l{TO|l for — 0, we 
find 

-^0"° = - ("pV +^'pV) k-^p:^pF (17) 

X™" = 2Ap o (upUp^ + VpVp') ~ 2Ap- 



with + Cp' = Lo and p -I- p' = k. 

The velocity of the collective mode — Vp/3 is small in the nonrelativistic 
system. Therefore, we expand the collective contribution in this parameter to 
obtain 

Mcon^ Ap i^^^l^^Py (18) 



cp 



The contribution of the bare vertex Mq'^'^'^ reproduces the matrix element 
((5J) derived by Flowers et al. ^ and Yakovlev et al. 0. However, the collective 
correction modifies this crucially. In the sum of the two contributions, the 
leading terms mutually cancel, yielding the matrix element 



Mo = M^"" + - o^F^ — 



which is at least ~ Vp times smaller than the bare result. 

The spatial component of the longitudinal (with respect to k) component 
of the matrix element can be obtained from Eq. H16|l . Since kF = {uj/k)To we 
have 

In the above, k = k/A; is a unit vector directed along the transferred momentum. 

Since the collective interaction modifies only the longitudinal part of the 
vertex, the transverse part of the matrix element can be evaluated directly from 
the bare vertex Hll|l . This yields 

M±~ [vpUp, - UpVp.) ^ ~ -^^1^^ (kp) Pi 

with p = p/p. 

The rate of the process is proportional to the square of the matrix element. 
This means that the vector current contribution to the neutrino energy losses is 
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Vp times smaller than estimated before. The corresponding neutrino emissivity 
in the vector channel can be evaluated with the aid of Fermi's golden rule: 



/f]3l. f]3l. 
^^c^Tr (l^Q M^M'-S (p + p' - k) <5 (cp + e^, 



One can simplify this equation by inserting J d'^K j''*-' (K — fci — /e2) = 1. Then, 
the phase-space integrals for neutrinos are readily done with the aid of Lenard's 
formula 

^^^(4) (X-fc,-fc,)Tr(Pn 

ArTT 

= — {K^K, - K'g^,) 6 {K^) 6 (lu) , 

where Q{x) is the Heaviside step function. 
For k <^ pp we obtain 

(2^) 

X ( (K^Mn^ - K^M^M^ ) S (2ep - u) 



The next integrations are trivial. We get 
592 



Qv = .oTo". ^ V^GlClpFM*T'y^ 



This is to be compared with Eq. We see that 

Qv 148 



2 
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Qfrs (i^°) 2835 



i.e. the neutrino radiation from ^Sq pairing in the nonrelativistic system {Vp <§; 1) 
is suppressed by several orders of magnitude with respect to the predictions of 
[3 and 13. 

Consider now the case when quasi-particles carry an electric charge. Includ- 
ing the long-range Coulomb interaction Vc (k) = e^/fc^ implies that the vertex 
part is multiplied by a string of closed loops, which represents the polarization 
of the surrounding medium. In this case, the new vertex F'^ can be found as 
the solution of the Dyson equation, according to the diagram of Fig. 1 or, 
analytically 

t''ip-K,p) = T^{p-K,p)^To{p-K,p)Vcik) 

d^p' 



Tr 



f^Qip' ~K)T^{p' -K,p')Q{p') (19) 
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Figure 1: Dyson equation for the vertex correction for charged quasi-particles. 
The shaded areas represent the modified effective vertex, and the wavy fine 
stands for the Coulomb interaction. 



This equation can be readily solved yielding 



[p IX, p) yp ix,p)yi l + Vc{k)W"{K)J 



with 



n°'^ {K)^t I Tr [t,Q {p' - K) [p' - K,p') Q {p')] . 

{2tt) 



In particular, for F*^ we arrive to^ 

The polarization function Ft"" (K) can be readily calculated with the help 
of F" given by Eq. (|14|l . By neglecting the small dependence of the energy gap 
on the transferred momentum k, we have 



n""(i^) = I 

• (2^) 
2uj 



J-0^Tr[r,Qip'-K)T,Qip') 



f^Tr [r,Q {p' - K) A [p') Q (p')] 
J {2n] 



- a'^p J (27r)'' 
Here, the quasi-particle propagator follows from Eq. (f7|: 



Qip) = - 2 (po+epT3 + fiA). (21) 

Po ~ 

We are interested in the regime defined by k < oj, uj > 2A ^ ep — ep_k — kVp- 
In this case we obtain, after some simplifications 

1 + X{K) 



•^The solution to the equation 1 + Vc (k) It"'' (K) = determines the new dispersion law 
Lu = u! (k) for the collective excitations, which represents plasma waves 1111 . 
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with 



(23) 



Since a ^ 1 and p ^ pp, the second integral in Eq. may be dropped. By 
neglecting also the small contributions from a?k^ ^ £pj'^^ we get 



xiK)^^-^ / ^^^^ 



with LOp = e^n/M* (n is the number of baryons per unit volume). This agrees 
with the plasma fr'equency for a fr'ee gas of charged particles. 

The energy exchange in the medium goes naturally as the temperature scale. 
Therefore, the energy transferred to the radiated neutrino-pair is w ^ T < 
Tc, while the plasma frequency ujp is typically much larger than the critical 
temperature for Cooper pairing. For instance, for a number density n of the 
order of the nuclear saturation density no — 0.17 fm and the effective mass of 
the baryon M* of the order of the bare nucleon mass, we obtain ujp ~ 10 MeV, 
while the critical temperature for baryon pairing is about 1 MeV or less. Under 
these conditions, we obtain 

f {p ^K,p)^ (p -K,p)^ [p - K,p) 

Thus, in superconductors, the vector current contribution to the neutrino radia- 
tion is suppressed additionally by a factor (T^ /^p) • this is the plasma screening 
effect. 

We have considered the problem of conservation of the vector weak current 
in the theory of neutrino-pair radiation from Cooper pairing in neutron stars. 
The correction to the vector weak vertex is calculated within the same order 
of approximation as the quasi-particle propagator is modified by the pairing 
interaction in the system. This correction restores the conservation of the vector 
weak current in the quasi-particle transition into a paired state. As a result, in 
the nonrelativistic baryon system, the matrix element of the vector current is 
Vp times smaller than previous estimations. This means that the vector weak 
current contribution to neutrino radiation caused by Cooper paring is Vp times 
smaller than it was thought before. The vector weak current contribution from 
pairing of charged baryons is suppressed additionally by a factor ~ (T^uip) 
due to plasma screening. The total suppression factor due to both the current 
conservation and the plasma effects is of the order 



iT!/u;iyV^< 10 
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Thus the neutrino energy losses due to singlet-state pairing of baryons can, 
in practice, be neglected in simulations of neutron star cooling. This makes 
unimportant the neutrino radiation from pairing of protons or hyperons. 
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